Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN WF$t style file v2.2) 



Non-Gaussian Foreground Residuals of the WMAP First 
Year Maps 



Joao Medeiros* and Carlo R. Contaldi 

Theoretical Physics, Imperial College, Prince Consort Road, London SW7 2BZ, UK 



5 February 2008 



ABSTRACT 

We inv estigate the effect of foreground residuals in the WMAP data l)Bennett et alJ 
2003a|) by adding foreground contamination to Gaussian ensembles of CMB signal 
and noise maps. We evaluate a set of non-Gaussian estimators on the contaminated 
ensembles to determine with what accuracy any residual in the data can be constrained 
using higher order statistics. We apply the estimators to the raw and cleaned Q, V, 
and W ba nd first year maps. Th e foreground subtraction method applied to clean 
the data in lBennett et al.l l)2003bf) appears to have induced a correlation between the 
power spectra and normalized bispectra of the maps which is absent in Gaussian 
simulations. It also appears to increase the correlation between the A£ = 1 inter-£ 
bispectrum of the cleaned maps and the foreground templates. In a number of cases 
the significance of the effect is above the 98% confidence level. 

Key words: cosmic microwave background - gaussianity tests. 



1 INTRODUCTION 

Over the past few years there has been a heightened inter- 
est in testing the statistical properties of Cosmic Microwave 
Background (CMB) data. The process has been accelerated 
by the release of the WMAP first year results. The WMAP 
data provide the first ever, full-sky maps which are signal 
dominated up to scales of a few degrees. Thus, for the first 
time we can test the Gaussianity and isotropy assumptions 
of the cosmological signal over large scales in the sample 
variance limit. 

Ever since the release of the COBE-DMR results 
(Benn ett et aljll99£li a consensus has been hard to reach 
on tests of non -Gaussianity with some studies report- 
ing null results (iKoeut et alJ Il996t IContaldi etafl l2000t 
ISandvik fc Magueiioll200*lF while others claimed detections 
of non-Gaussian features dFerre ira et al ] ll99St iMaeueiic 
20001: iNovikov. Feldman fc Shandarinl Il998i: iPando et al 
19981) 



With the release of the WMAP first year results 
a limit on the non-Gaussianity of primordial perturba- 
tions in the form of an estimate of the non-linear factor 



fm, was obtained bv iKomatsu et alJ J2003I 


) . However a 


number of authors (IBielewicz et alJ 12004: 


Erikscn et alJ 


2004al: IColes et alJ 12004]: iParkl 


120041: ICopi et all 


2004; 


de Oliveira-Costa et all boo4 


Land & Maeueiia 


20051; 


Efstathioull2003l iRoukema et alJ 


20041: lEriksen et all 


2005; 


Jaffe et al. 200lf) have also reported analysis 


of the maps 
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that suggest violations of the Gaussian or isotropic nature 
of the signal. 

One of problems with testing Gaussianity is that one 
can devise a plethora of tests to probe the infinite degrees 
of non-Gaussianity, therefore different tests represent dif- 
ferent perspectives on the statistical patterns of the signal. 
For WMAP there are already a number of detections of 
so called anomalies, most pointing to different unexpected 
features in the mi c rowave sky. The most documented case 
JPeiris et ai]l2003t lEfstathioJ l2003t ISlosar fc Seliakl 120041 : 
Ide Oliveira-Costa et all 120041: lEriksen et alJ l2004af) is the 
low amplitude of the quadrupole and octupole in compar- 
ison to the inflationary prediction, something we can cat- 
egorize as amplitude anomalies. Although it is simple to 
design inflationary spectra with sharp features which repro- 
duce, more or l ess clo se ly, the ampli t ude a n omaly (see e.g. 
IContaldi etafl <2003ft : iBridle et al] (l2003ft ; ISalopek et all 
l^89Tv Tthese invariably suffer fine tuning problems. Another 
approach is to relate the anomaly to the breakdown of sta- 
tistical isotropy or Gaussianity. 

Other reported features relate to the correlation of 
phases in the multipole coefficients which are an indication 
of non-Gaussianity. These can be dubbed ph ase anomalies. 
One e xample is the hemisphere asymmetries (Erikscn et al. 
l2004al) : the northern ecliptic hemisphere is practically flat 
while the southern hemisphere displays relatively high fluc- 
tuations in t he power spectrum. Othe r functions, such as the 
bispectrum jLand fc Magueij ol l2005h and n-point correla- 
tion functions 7ETiksen^t^!Tl l2005 j also show related asym- 
metries. Furthermore, there is the anomalous morphology 
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of several multipoles, in particular, the striking planarity 
of the quadrupole and octupole a nd the strong alignment 
betw een their preferred directions l|de Oliveira-Costa et alJ 
12004) . Overall, there is a strong motivation to continue prob- 
ing the statistical properties of the data and find possible 
sources for these signals, be it instrumental, astrophysical 
or cosmological. 

The first test to have provided indications of possible 
non-Gaussian featur es in the CMB data was reported by 
iFerreira et, alJ (|f99St) and IMagueiiol l)2000Fl using a bispec- 
trum estimator, the fourier analog of the three point func- 
tion. Both those detections were later found to be caused 
by systemat ic effects rathe r than by cosmological source as 
reported by iBandav et alJ J200(J) . For the ca se of the bis- 
pectrum signal detected bv IMagueiiol feOOOl) . which used 
an estimator tuned to detect correlations between neigh- 
bouring angular scales, finding the source of the signal had 
to wait for the releas e of the high precision WMAP data 
^Bennett et all2003aT) which was able to provide a compara- 
tive test of the cosmological signal. The WMAP data did not 
reproduce COB E's result and systematic err ors were found 
to be the cause l|Magueiio fc Medeirosl 12004 ). The WMAP 
da ta was later analysed wit h the bispectrum in more detail 
bv lLand fc Magueiiol J2005T) . In that paper, the bispectrum 
of the clean, coadded maps was analysed and a connection 
between the hemisphere asymmetries in the 3-point correla- 
tion function and the bispectrum was established, although 
the full sky as a whole was found to be consistent with Gaus- 
sianity. 

In this paper, we study the effect that foreground con- 
taminations have on bispectrum estimators. In section|5|we 
define a set of bispectrum estimators with set I configura- 
tions. In section 13 we describe the template dust, free-free 
and synchrotron maps used to characterize the effect on the 
bispectrum. In section|l]we determine the distribution of the 
estimators in the presence of residual foregrounds with dif- 
ferent amplitudes and discuss the application of this method 
to detect residuals in the data by introducing a number of 
statistical and correlation measures. In section[2]we discuss 
the application of the the statistical tools developed in the 
previous sections to the raw and cleaned WMAP first year 
maps. We conclude with a discussion of our method and 
results in section^] 



2 THE ANGULAR BISPECTRUM 

We now introduce the angular bispectrum estimator 
iFerreira et alj l"l998). The bispectrum is related to the third 
order moment of the spherical harmonic coefficients at m of a 
temperature fluctuation map AT(n)/T. The coefficients de- 
scribe the usual expansion of the map over the set of spher- 
ical harmonics Y( m (h) as 



AT 

— (n) - ^ 

Ira 



aimYe m (n). 



(1) 



Given a map, either in pixel space or harmonic space, 
and assuming statistical isotropy, one can construct a set hi- 
erarchy of rotationally invariant statistical quantities char- 
acterizing the pattern of fluctuations in the maps. These are 
the n-point correlation functions in the temperature fluctu- 



ations {=p-(m)=^-(n)...=p-(p)} or in the spherical harmonic 
coefficients, {a£ imi ai 2m2 ■■■ae n m n }- 

The unique quadratic invariant is the angular power 
spectrum defined as {ai imi at 2 m 2 ) = S ei i 2 S mim2 C e , whose 
estimator can be written as Ce = ojqrr X/rr l a *™l 2 - This 



gives a meas ure of the ove r all int ensity for each multipole 
£. Following IFerreira et alJ il99STl . the most general cubic 
invariant defines the angle averaged bispectrum, 
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where the (. . .) is the Wigner 3J symbol. Parity invariance 
of the spherical harmonic functions dictates that the bispec- 
trum be non-zero only for multipole combinations where the 
sum £1+ £2 + tz is even. An unbiased estimator (for the full 
sky) can be evaluated as 
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with the normalization factor defined as 







(4) 



(2£i + l) (2^ + 1) (24 + 1) 
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The bispectrum can be related to the three-point correlation 
functions of the map just as the power spectrum Ct can be 
related to the correlation function C(0) through the well 
known expression 



C(0) 



-y 



(2£+l)C^(cos0). 



(•») 



For example, the pseudo-collapsed, three-point correlation 
function, C (3) (<9) = (^(n) 2 ^(m)}, is related to our defi- 
nition of the bispectrum Bi 1 e 2 e 3 as 



(6) 



where n • rh = cos 6. 

It is important to use both tools, the bispectrum and 
the three-point correlation function, to probe the sky maps 
as they have the capacity to highlight different features of 
the data. In principle, harmonic space based methods are 
preferred for the study of primordial fluctuations whereas 
real space methods are more sensitive to systematics and 
foregrounds, which are strongly localized in real space. In 
addition, the three-point correlation function is intrinsically 
very sensitive to the \aw-£ modes, whereas the bispectrum 
can pick up different degrees of freedom with respect to the 
different mode correlations we want to probe 

For the choice £1 = £2 = £3 = £ we can de fine the 
single-^ bispectrum Be — Bin iFerreira et aljll998f) . which 
probes correlations between different m's. Other bispectrum 
components are sensitive to correlations between different 
scales £. This can be extended to study correlations from 
different angular scales. The simplest of these is the A£ — 1 
inter-£ bispec trum between ne ighbouring multipoles defined 
as Bt-itn+i iMagueiiol '20001 It is convenient to consider 
estimators normalized by their expected Gaussian variance 
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Ce 1 Ci 2 Ce 3 which have been shown to be more optimal and 
more Gaussian distributed than the unnormalized estima- 
tors, and are not sensitive to the overall power in the maps. 
Here we will introduce the Ie,Je, and K e bispectra defined 



B, 



and 



K 



(ft) 3 / 2 



Bp 



T } — 



1 et+i 



(7) 



(8) 



(Cl-2ClCt+2) 1/2 

where have extended the formalism to a separation Al = 2 
to probe signals with both odd and even parity in the inter-£ 
correlations. 



3 FOREGROUND TEMPLATES 

The standard method of foreground removal used by cos- 
mologists makes use of a set of template maps for each of 
the dominant sources of foreground contamination in the 
CMB frequency maps. These are maps obtained from in- 
dependent astronomical full-sky observations at frequencies 
where the respective mechanisms of emission are supposed to 
be do minant. These templates are the H a map llFinkbeinerl 
|2003j), for the free- fr ee emission, the 408 MHz Haslam map 
jHaslam et"ai]|l98ll) . f or the synchrotron e mission, and the 
FDS 94 GHz dust map JSchleeel et al.ll998l) . These are then 
subtracted from the WMAP data with coupling coefficients 
determined by cross correlating with the observed maps in 
the Q (41 GHz), V (61 GHz), and W (94 GHz) bands. Nev- 
ertheless the templates are a poor approximation of the of 
the real sky near the galactic plane, so a Kp2 mask must 
still be used in the a naly sis. The method is de scribed in 
iBennett etafl J2003bl) and lKomatsu et alJ J2002h : 
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where r\ is a correction factor due to reddening in the free- 
free template and vq = 40.7 GHz, vq = 60.8 GHz and 
v w = 93.5 GHz . The values in front of the left bracket 
convert the detector's temperature to thermodynamic tem- 
perature. It is considered that this is a sufficiently good 
method to remove the foregrounds outside the Kp2 plane 
since it matches the correct amplitudes quite well, however 
the usual doubts remain, especially in the light of the align- 
ment/low multipoles controversies. Another point one can 
make is that whereas this may be a satisfactory technique 
to correct the foregrounds at the power spectrum level, its 



effect on higher order statistics is unknown and may actually 
induce unexpected correlations. 



4 THE EFFECT OF FOREGROUNDS ON THE 
BISPECTRUM 

We have generated a set of 3000 Gaussian, CMB simu- 
lations of t he WMAP first y ear Q, V, and W maps in 
HEALPix 1 iGorski et al.ll2005i) format with a resolution pa- 
rameter iV s ido = 512. Each simulation is smoothed with the 
Q, V and W frequency channel beams and channel spe- 
cific noise is added. We adopted the WMAP best-fit ACDM 
with running index power spectrum 2 to generate the ai m 
coefficients of the maps. The Kp2 galactic mask is imposed 
on each map. The masked maps are then decomposed into 
spherical harmonic coefficients at_ m using the Anafast rou- 
tine. We then calculate the four spectra; namely the the 
power spectrum Ci, single-^ bispectrum if, A£ = 1 inter-^ 
bispectrum jf and A£ — 2 inter-^ bispectrum Kf as de- 
scribed in section [5] 

We then add channel-specific foregrounds outside the 
Kp2 zone to the same set of Gaussian simulations with am- 
plitudes set as in Eqn. @. The addition of the foreground 
is scaled linearly by a factor a as 



T{Q,v,w} = Tcmb + aT 



{Q,V,W}, 



(10) 



which we use to check the sensitivity of the bispectra to 
the foregrounds (typically a = 1.0 or a = 0.5). The power 
spectrum and bispectra are then calculated for the set of 
contaminated maps. 

In Fig.^we show the mean angular spectra of the simu- 
lations obtained by averaging over the ensembles. We show 
the mean spectra for the Gaussian (solid, black) and the 
contaminated simulations for a — 0.5 (short-dashed, red) 
and a — 1.0 (long-dashed, blue). The shaded area shows the 
variance of the three bispectra obtained directly from the 
Gaussian simulations. 

We see that even for the fully contaminated set of maps 
(q = 1.0) the average signal is not significantly larger than 
the expected Gaussian variance indicating that a detection 
would require averaging over a large number of modes. How- 
ever we see some important distinguishing features in the 
signal in that it is sensitive to the parity of the multipole, 
being suppressed for odd I. This is due to the approximate 
symmetry of the foreground emission about the galactic 
plane which means that most of the signal will be in even I 
modes since these have the same symmetry. This effect can 
be seen in all the spectra but most significant is the sup- 
pression of the odd inter-i 1 bispectrum jf with respect to 
the even inter-£ bispectra if and K\. 

Another obvious feature of the even parity nature of 
the signal is the correlation between the spectra. In partic- 
ular the absolute values of the if and the Kf are correlated 
with the structure visible in the fully contaminated power 
spectrum. 

Overall the Kf is the most sensitive statistic with the 
largest amplitude with respect to the Gaussian variance al- 
though still quite small even at 50% contamination. We now 



1 http://healpix.jpl.nasa.gov 

2 http://lambda.gsfc.nasa.gov 
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Figure 1. The average functions for the angular spectra of the simulations. Black (solid) is Gaussian, red (short-dashed) is for the 
contaminated simulations with a = 0.5 and blue (long-dashed) is for the contaminated simulations with a = 1.0. The top panel shows 
the power spectrum Ci , second panel shows the single- £ bispectrum If, the third panel the A£ = 1 inter-^ bispectrum Jf and the bottom 
panel shows the A£ = 2 inter-£ bispectrum Kf . The shaded regions represent the Gaussian variance measured directly from the ensemble 
of simulations. For the Gaussian simulations, the average power spectrum is just the input ACDM power spectrum and the average 
bispectra is effectively zero. On the other hand, the average angular spectra of the contaminated simulations have an emerging pattern 
of intermittency in both second- and third-order statistics. This intermittent pattern comes about due to the even parity of galactic 
foregrounds, ie, even modes are enhanced relatively to the odd modes. This can be seen in the significant increase of power in the even 
modes of the power spectrum. In terms of the bispectrum, we see that the A£ = and the A£ = 2 inter-4? components will be more 
significantly enhanced than the A£ = 1 inter-f bispectrum because the latter includes correlations between even and odd modes. 



describe a number of statistical estimators we use to test the 
detectability of the template matched foregrounds in the Q, 
V, and W channel maps. 

4.1 Chi-Squared Test 

Having seen how foregrounds affect the angular statistics of 
CMB maps, we can now devise specific tests to probe these 
properties on the bispectrum and test their sensitivity. The 
standard way to use the bispectrum as a tes t of general non- 
Gaussianity is to use a reduced y 2 statistic IIMagueiioll2000l: 
iMaeueiio fc Medeirosl20ol lLand fc MagueiicfeOOrJ) . This is 
defined as 
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(X t - {X,)f 



(11) 



where X( is a given bispectrum statistic, (Xt) is its mean 
value computed over the Monte Carlo ensembles, and erf is 
the variance for each angular scale. The \ 2 test is a mea- 
sure of the deviation of the observed data from the expected 
mean, weighted by the Gaussian variance of the estimator. 



Foregrounds increase the amplitude of the bispectra 
foregrounds, but as shown in Fig. we can see that only 
Kf seems to stand of chance of significant detections since 
the average amplitude of the signal is comparable to the 
variance, unlike the other components of the bispectrum. 

The detectability of the template matched signals using 
any of the bispectra can be tested by comparing the dis- 
tribution of the x 2 values obtained from the contaminated 
simulations with that obtained from Gaussian simulations. 
We compute the \ 2 values for the contaminated maps us- 
ing the mean and the variance obtained from the Gaussian 
simulations, ie, the expected Gaussian functions. 

We compare the distribution of the \ 2 values for the 
Gaussian simulations against the distribution obtained for 
the simulations with contamination (a = 1.0). We concen- 
trate on the Q band since it is the most contaminated fre- 
quency. The histograms of the \ 2 are shown in the left col- 
umn of Fig. [5] For the if and jf spectra the histograms 
overlap completely. This means that the probability of find- 
ing contaminated simulations with a high x 2 is the same as 
for the Gaussian simulations indicating that the \ 2 test is 
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insensitive to the presence of foreground contaminations at 
this level. However the K\ spectrum tells a different story. 
There is a significant shift between the two distributions 
which implies that this component of the bispectrum has 
more sensitivity to foregrounds. 

The sensitivity can be quantified in terms of the frac- 
tion of the contaminated simulations (with a = 1.) with a 
\ 2 larger value than 95.45 % (i.e. 2 a) of the Gaussian sim- 
ulations (with a = 0). The sensitivity for the if and Jf is 
low (i 0.05), whereas for K\ the fraction increases to 0.355. 



4.2 Template Correlation Test 

A template matched statistic can be defined by correlating 
the observed bispectra in the data with those of the fore- 
ground templates. This is more sensitive to the structure 
in the template signal as opposed to the \ 2 test introduced 
above. We define a cross correlation statistic p as 



XiX £ 



Efmax y2 V' ma 



Xf 



1/2 



(12) 



where Xt are the bispectra obtained from the data and the 
Xf are those obtained from the foreground templates. 

In the middle column of Fig.|5]we display the histograms 
for the p values for the Gaussian simulations against the dis- 
tribution obtained for the contaminated (a = 1) simulations 
of the Q band maps. The sensitivity has improved over the 
X 2 test, with the histograms of the input and output data 
sets being clearly shifted, meaning that there is a higher 
probability of detection of foregrounds using this method. 
Again the effect is stronger in the K\. This result simply 
quantifies the statement that a matched template search for 
a contamination signal is more sensitive than a 'blind' statis- 
tic such as the \ 2 test. The values for the sensitivity of the 
test are given in tableQfor all three WMAP bands. 



4.3 Power Spectrum and Bispectra 
Cross-Correlation Test 

For a Gaussian field, the normalized bispectrum is statis- 
ticall y uncorrelated with the power spectrum jMagueiid 
Il995l) . However, foreground residuals in the map induce non- 
Gaussian correlations which in turn will induce correlations 
between the normalized bispectra and the power spectrum 
of the maps. This can provide another specific signature that 
one can use to detect the presence of foreground contami- 
nation. 

For Gaussian simulations, the average power spectrum 
is just the input ACDM power spectrum and the bispec- 
trum is effectively zero. On the other hand, the average 
angular spectra of the contaminated simulations have an 
emerging pattern of intermittency in both first- and second- 
order statistics. Correlations between the power spectrum 
and the bispectra therefore come about due to the even par- 
ity induced by the characteristic galactic foregrounds. This 
means that the even modes of the power spectrum will be 
correlated with the even modes of the bispectra, whereas 
odd modes will remain uncorrelated. In order to test this 
effect on the maps, we introduce the R correlation statistic 
defined as 



Table 1. Sensitivity of the p and R tests in terms of the fraction 
of the contaminated simulations (a = 1.0) with a larger value 
than 95.45 %, ie 2 <r, of the Gaussian simulations (a = 1.0). We 
present values for the Q,V and W frequency channels and for 
the I?, J? and Kf. Note that p(JQ), where X( is a given bis- 
pectrum component stands for the correlation between Xg (data) 
withJQ (template), whereas R(JQ), represent the correlation of 
that specific bispectrum component with the respective power 
spectrum of the map. The values in the table quantify what can 
be seen in the histograms in figure [2] Applying the tests for the 
Kf component provides better sensitivity to the foregrounds. Be- 
tween the two tests p seems to provide a marginally better sensi- 
tivity. Also, the Q channel, being the most foreground contami- 
nated yields the higher chances of detection. 
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Table 2. Results for the WMAP data for p and R. The results 
are shown as the fraction of Gaussian simulations below the level 
obseverved in the data. We have highlighted values with greater 
than 98% in the foreground cleaned maps. 
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where Ci is the observed power spectrum. We have chosen 
^max = 30 and £ m i n = 4 as we are interested in the large 
angular scales where the effects of foreground contamination 
will dominate. We use the absolute value of the bispectrum 
in order to avoid the discrimination between negative and 
positive correlations which would affect our sum. We are 
only interested in the discrimination between the existence 
of absolute correlations against null correlations between the 
Ce and the bispectra Xe . 

Again we test the sensitivity of this method by com- 
puting a distribution of R for Gaussian ensembles against 
the contaminated ensembles. We make sure that for Gaus- 
sian ensembles we use the correlation of Cf +F with Xf and 
for the contaminated ensemble the correlation of C^ +F with 
Xf +F where S stands for the Gaussian CMB signal and S+F 
indicates contaminated ensembles. This allows us to cancel 
the effect of the increase of power due to foregrounds in the 
correlation of the two statistics between the two tests. The 
results for the contaminated ensemble, a = 1, are plotted in 
the right column of Fig.|2|and are summarized in tablcQfor 
all three bands. 
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Figure 2. Distributions of values obtained for the three different tests applied to the Q channel to detect the presence of foregrounds 
(x 2 , p and R), for Gaussian simulations (black) and contaminated (a = 1.0) simulations (grey). The level of sensitivity of a given method 
can be determined in terms of the shift between the histograms for the Gaussian and the contaminated case. We see that the p K and 
R^ statistics are the most sensitive channels to probe the existence of foregrounds. 



5 APPLICATION TO THE WMAP DATA 

We have applied the stat istical tools described above to the 
WMAP first year data jBennett et alJl2003aT) . We consid- 
ered both the raw and cleaned maps of the Q, V, and W 
channels using the Kp2 exclusion mask. We summarise the 
results in table|5|showing the separate confidence limits from 
each channel for both the raw and cleaned maps. 

For the raw maps we find that only the Q channel R x 
result is above the 95% threshold while for the Q channel p 
statistic, all confidence levels are above the 90% level with 
the p 1 above the 95% level. This is consistent with there 
being a component most correlated to the foreground tem- 
plates at the lowest frequencies and with significant correla- 
tions between the A£ = 2 inter-^ bispectrum and the power 
spectrum. Since the raw maps do not have any foreground 
subtracted from them this is not a surprise although the con- 
fidence level suggests that the correlations are larger than 
what was found for the expected amplitude (a — 1) of the 
foregrounds. 

For all if and Kf statistics the cleaned map results 
show confidence levels below the 95% level and indeed show 
an overall reduction in the significance of the correlations, 
indicating that the cleaning has removed a component cor- 
related to the foreground templates, as one would expect. 
However for the jf statistics, which should in principle be 
the least sensitive to the foregrounds considered, we see that 
the confidence levels have all increased. Indeed all three 



channels now have correlations significant above the 95% 
level in the p statistic with the W channel also having 
a > 95% confidence level. The cleaning algorithm appears to 
have introduced significant correlations with the foreground 
templates in the A£ — 1 inter-£ bispectra and significant cor- 
relations between the At — 1 inter-^ bispectrum and power 
spectrum of the W channel which is indicative of a non- 
Gaussian component. 

In figure[3]we show the bispectra for each cleaned chan- 
nel map and compare to the bispectra of the foreground 
template (a = 1) for each channel. This shows the nature 
of the result above. For both the if and Kf the cleaned 
map bispectra are anti-correlated with the foreground tem- 
plates. In addition the the Kf for all channels are heavily 
suppressed in the cleaned maps for multipoles £ < 20 com- 
pared to the expected Gaussian variance shown in figure 
The Jf gives the only bispectra that are correlated with the 
those of the templates. 

Figure^Jshows the break down of the R J result into in- 
dividual multipole contributions for each of the three bands. 
In particular it is interesting to note how the W band R J 
result shown in tabled is dominated by an outlier at £ — 26. 



6 DISCUSSION 

At first sight our results appear contradictory. We have stud- 
ied the effect of foreground contamination on the maps and 
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Figure 3. The bispectra for each cleaned channel map (solid line) against the bispectra of the foreground template (a = 1.0) (dashed 
lines). The three frequency channels are shown as Q (red), V (black) and W (blue). Both the cleaned map if and Kf bispectra are 
anti-correlated with the foreground templates. In addition the the Kf for all channels are heavily suppressed in the cleaned maps for 
multipoles i < 20 compared to the expected Gaussian variance shown in figure Q The jf gives the only bispectra that are correlated 
with the those of the templates. 



concluded that foregrounds mainly affect the if and Kf 
components of the bispectrum due to its parity. By com- 
paring the results for the raw and the foreground-cleaned 
maps, we are able to verify that the amplitude of if and Kf 
reduces as expected after foreground subtraction. 

On the other hand, as shown in table [5] the correla- 
tions induced in the jf appear to be close to inconsistent 
to a Gaussian hypothesis with the correlation with the fore- 
ground templates at a significance above the 3<r level for the 
Q-band, cleaned map. It is also of interest to note that the 
cleaned maps do worse in all bands for the p measure. 

This is not what we naively expected since the fore- 
grounds considered here have the wrong parity and their 
jf signal is heavily suppressed. However the cleaning pro- 
cedure used by the WMAP team does appear to increase 
the correlations p of jf bispectrum to the input maps and 
its correlation R with the power spectrum. Recall that we 
expect the normalized bispectra to be independent of the 
power spectrum only in the Gaussian case. 

The possibility of the foregrounds being more complex 
than accounted for in this type of treatment is to be con- 
sidered carefully as this work has shown. The results shown 
here would suggest that the procedure used to go from the 
raw to cleaned WMAP maps is under or over subtracting 
a component with £ ± 1 parity in the bispectrum. This is 
probably not an indication that the procedure is faulty but 



rather that the templates used are not accurate enough to 
subtract the foregrounds. One source of inaccuracy is the 
simple scaling of the templates with respect to frequency. 
The clean ed maps are ob t ained assuming uniform spectral 
index and iBennett et alJ <l2003bl) acknowledge that this is 
a bad approximation particularly for the 408 MHz Haslam 
(synchrotron) template. This is seen when producing the 
Internal Linear Combination (ILC) map which accounts for 
variation of the spectral index of the various component. Un- 
fortunately ILC maps cannot be used in quantitative studies 
as their noise attributes are complicated by the fitting pro- 
cedure and one cannot simulate them accurately. 

Future WMAP ILC maps or eq uivalent ones obtained 
by 'blind' foregroun d subtraction llTegmark et alJ 120031 : 
lEriksen et al.ll2004bl) may be better suited for this kind of 
analysis once their statistical properties are well determined. 
It is expected that the impending second release of WMAP 
data will allow more accurate foreground analysis and the 
statistical tools outlined in this work will be useful in deter- 
mining the success of foreground subtraction. 

It may be worthwile to include information of the higher 
order statistics when carrying out the foreground subtrac- 
tion itself, for example by extending the ILC method to 
minimise higher order map quantities such as the skewness 
and kurtosis of the maps. 
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Figure 4. R J as a function of angular scale, £. The values displayed correspond to the foreground cleaned Q (red squares), V (black 
triangles) and W (blue circles) frequency channels. The values are offset by i = 0.25 and I = 0.5 for the Q and W bands respectively. 
The 3cr detection in the W foreground-cleaned channel is dominated mainly by the I = 26 mode. The error bars are computed from 3000 
Gaussian simulations assuming the specific channel noise and beam. 
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